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ABSTRACT 

Following the seminal result of An & Evans, known as the central density slope- 
anisotropy theorem, successive investigations unexpectedly revealed that the density 
slope-anisotropy inequality holds not only at the center, but at all radii in a very large 
class of spherical systems whenever the phase-space distribution function is positive. 
In this paper we derive a criterion that holds for all spherical systems in which the 
augmented density is a separable function of radius and potential: this new finding 
allows to unify all the previous results in a very elegant way, and opens the way for more 
general investigations. As a first application, we prove that the global density slope- 
anisotropy inequality is also satisfied by all the explored additional families of multi- 
component stellar systems. The present results, and the absence of known counter- 
examples, lead us to conjecture that the global density slope-anisotropy inequality 
could actually be a universal property of spherical systems with positive distribution 
function. 
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1 INTRODUCTION 

In the study of stellar systems based on the "p-to-/" ap- 
proach (where p is the density distribution and / is the as- 
sociated phase-space distribution function, hereafter DF), 
p is given, and specific assumptions on the internal dynam- 
ics of the model are made (e.g. see Bertin 2000; Binney & 
Tremaine 2008, hereafter BT08). In some special cases inver- 
sion formulae exist and the DF can be recovered in integral 
form or as series expansion (see, e.g., Fricke 1952; Lynden- 
Bell 1962; Osipkov 1979; Merritt 1985; Dejonghe 1986, 1987; 
Cuddeford 1991; Hunter & Qian 1993; Ciotti & Bertin 2005). 
Once the DF of the system is known, a non-negativity check 
should be performed, and in case of failure the model must 
be discarded as unphysical, even if it provides a satisfac- 
tory description of data. Indeed, a minimal but essential 
requirement to be met by the DF (of each component) of 
a stellar dynamical model is positivity over the accessible 
phase-space. This requirement, the so-called phase-space 
consistency, is much weaker than the model stability, but 
it is stronger than the fact that the Jeans equations have a 
physically acceptable solution. However, the difficulties in- 
herent in the operation of recovering analytically the DF 
prevent in general a simple consistency analysis. 
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Fortunately, in special circumstances phase-space con- 
sistency can be investigated without an explicit recovery 
of the DF. For example, analytical necessary and sufficient 
conditions for consistency of spherically symmetric multi- 
component systems with Osipkov-Merritt anisotropy (Os- 
ipkov 1979, Merritt 1985, hereafter OM) were derived in 
Ciotti & Pellegrini (1992, hereafter CP92; see also Tremaine 
et al. 1994) and applied in several investigations (e.g., Ciotti 
1996, 1999; Ciotti & Lanzoni 1997; Ciotti & Morganti 2009, 
hereafter CM09; Ciotti, Morganti & de Zeeuw 2009). More- 
over, in Ciotti & Morganti (2010, hereafter CMlOa) ana- 
lytical necessary and sufficient consistency criteria for the 
family of spherically symmetric, multi-component, general- 
ized Cuddeford systems (which contain as very special cases 
constant anisotropy and OM systems) have been derived. 

Another necessary condition for consistency, the focus 
of this paper, is the "central cusp-anisotropy theorem" (An 
& Evans 2006, hereafter AE06; see also equation [28] in de 
Bruijne et al. 1996), an inequality relating the values of the 
central density slope 7 and of the central anisotropy pa- 
rameter P of consistent spherical systems, namely 7 > 2/3 
(see Section 2). This condition was derived for constant 
anisotropy systems, and then generalized asymptotically 
to the central regions of spherical systems with arbitrary 
anisotropy distribution. In AE06 it was also shown that 
the density slope-anisotropy inequality actually holds rig- 
orously at every radius in constant anisotropy systems, and 
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not only at their center. We will refer to this case, i.e. when 
y(r) > 2/3(7-) Vr, as the Global Density-Slope Anisotropy 
Inequality (hereafter GDSAI). 

Surprisingly, in CM09 we showed that the CP92 neces- 
sary condition for model consistency is nothing else than the 
GDSAI in disguise. Thus, not only in constant anisotropy 
systems but also in each component of multi-component 
OM systems the GDSAI holds. Prompted by this curious 
result, in CMlOa we introduced the larger family of multi- 
component generalized Cuddeford systems, we studied their 
phase-space consistency, and we finally proved that the GD- 
SAI is again a necessary condition for phase-space consis- 
tency of each density component. 

The results of CM09 and CMlOa revealed an unex- 
pected generality of the GDSAI and, in absence of known 
counter-examples (see also the discussion in CMlOa), it is 
natural to ask whether the GDSAI is even more general, 
i.e. it is necessarily obeyed by all spherically symmetric, 
two-integrals systems with positive DF. If such conjecture 
proved true, it would be remarkable not only from a theo- 
retical point of view, but also for applications. In fact, as 
it would hold separately for each density component of a 
stellar system, the value of the anisotropy parameter of the 
stellar component would be controlled at each radius by the 
local stellar density slope itself, independently of the dark 
matter halo. This constraint could be then used to reduce 
the impact of mass-anisotropy degeneracy in observational 
work^J. Motivated by the above arguments, we searched for 
a proof of the general validity of the GDSAI. While a proof 
is still missing, some relevant advance has been made: in 
particular, we obtained a new criterion that allows us not 
only to prove in a very elegant and unified way all our pre- 
vious results, but also to demonstrate that new families of 
models do necessarily obey the GDSAI when consistent. 

The paper is organized as follows. In Section 2 a general 
criterion linking phase-space consistency to the GDSAI for 
systems whose augmented density is a separable function of 
radius and potential is derived, which unifies all the results 
obtained so far on the GDSAI, and opens the way to the 
investigation of an even wider class of models. In Section 3 
the new criterion is used to prove in a new way that the GD- 
SAI is obeyed by generalized Cuddeford models, but also by 
some well-known stellar systems not belonging to the fam- 
ily of generalized Cuddeford models, thus further extend- 
ing the validity of the GDSAI as a necessary condition for 
phase-space consistency. Finally, the main conclusions are 
summarized in Section4. 



2 A GENERAL CRITERION 

In Ciotti & Morganti (2010, hereafter CMlOb) we showed 
analytically, by direct computation, that two well-known 
anisotropic models not belonging to the generalized Cudde- 
ford family, namely the Dejonghe (1987) anisotropic Plum- 
mer model, and the Baes & Dejonghe (2002) anisotropic 
Hernquist model, indeed obey the GDSAI whenever their 

1 Clearly, as by definition /3 < 1, such limitation would be useful 
only in the regions where the density slope is < 2, while in the 
galaxy outskirts the inequality is not helpful, as 7 > 2 for mass 
convergence. 



DF is positive. On one hand this result proves that the global 
inequality is not a specific property of generalized Cudde- 
ford models, due perhaps to some special dependence of their 
DFs on energy and angular momentum. On the other hand, 
it reinforces the conjecture that the GDSAI could be a very 
general (if not a universal) property of consistent spherical 
models. In this Section we provide a new hint to the latter 
hypothesis, as we show, with the aid of a new criterion, how 
the GDSAI is rigorously valid in a very large class of consis- 
tent systems, containing not only the multi-component gen- 
eralized Cuddeford systems and the two models discussed in 
CMlOb, but whole new families of models. 



2.1 General relations and the case of factorized 
systems 

We consider stationary, non-rotating, spherically symmet- 
ric systems with a two-integrals phase-space distribution 
function / = /(£, J), where £ = — v 2 /2 is the relative 
energy per unit mass, = —&t is the relative total po- 
tential, and J is the angular momentum modulus per unit 
mass. In general may contain also the contribution of an 
"external" potential (for instance the one corresponding to 
a dark matter halo). 

It is easy to show (e.g. see BT08, Ciotti 2000) that the 
density distribution p, the radial velocity dispersion <r r , and 
the tangential velocity dispersion er t are related to the DF 



P = % I d£ 



Jm f(£, J) J 



dJ, 
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d£ / f(£,J)AJdJ =p T , 
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f(£,J)J 3 



dJ = p t , 



(1) 
(2) 
(3) 



i) 



where A = V /2(* T - £) ~ J 2 /r 2 and J m = r v /2(* T - £ )■ 
With p T and p t we indicate the radial and tangential pres- 
sure, respectively. 

Given the identities above, and given the definition of 
the anisotropy parameter f3, easy algebra proves the remark- 
able identities 

dp. 



p(r,* T ) 
0(r,*T) 



1 



Pt 

2 Pl 



1 d\np t 

2 dlnr 



(4) 
(5) 



(e.g., see Spies & Nelson 1974, Dejonghe 1986, Cavenago 
1987, Dejonghe & Merritt 1992, Bertin et al. 1994, Bertin 
2000, Baes & van Hese 2007). Note that these relations hold 
independently of the specific radial dependence of $t, and 
that the radial trends of p and /3 are known only after the 
total potential >Pt is given. The existence of the general 
relations (Q and (5) is of the greatest importance for the 
present study, because it shows clearly that p (and so 7) 
and P are somewhat linked by the function p T . This link 
could open the way to a general proof of the GDSAI. 

At the present stage, we do not attempt a general proof, 
but we focus our attention on a special case, that of spherical 
systems with a factorized augmented density. In practice, 
we consider special spherical systems in which the radial 
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pressure p r is a factorized function of the radius and of the 
total potential. Therefore, from equation Q, it follows that 
also the density distribution can be written as 



P (r,* T )sA(r)B($, 



(6) 



of course, while the function B in the expression above is 
the derivative of the potential dependent factor in the fac- 
torized expression of p r , the radial function A is the same. 
We note that the multi-component generalized Cuddeford 
models (and therefore also the constant anisotropy, OM, and 
Cuddeford models, see CMlOa), as well as the two models 
discussed in CMlOb, belong to such family of systems. From 
equation (5) it follows that, independently of the specific ra- 
dial dependence of fx, 

2 » - -Sf- m 

Now, for assigned ^t{t), equation §4$ shows that the 
logarithmic slope of the density profile (|6} can be written as 

dlnp dlnB dln* T 



7(r) 



din r 



din * T din r 



+ 2/3(r); 



(8) 



since from Newton theorem Vt^r) is a monotonically de- 
creasing function of radius, identity ([8j proves the following 
Criterion 1: in all spherical systems whose density dis- 
tribution is a separable function of radius and total poten- 
tial, p = A(r)_B(# T ), the global inequality -y(r) > 2j3{r) Vr 
holds B('I't) is a monotonically increasing function of 

Therefore, if one is able to show that in all factorized 
consistent systems the B function is necessarily monotonic, 
then the GDSAI will hold in all these systems. We were not 
able to prove or disprove this possibility in general, however 
in Section 3 we present interesting results along this line. 

Before moving to discuss the new results, we note the 
following curious fact: the GDSAI can also be expressed as 
a condition on the radial velocity dispersion of the model. In 
fact, for a two-integrals spherical system, the relevant Jeans 
equation can be written as 



dpa? i 2/3 pa? 
dr 



+ 



d* T 
dr 



(9) 



(e.g., BT08). Introducing the logarithmic density slope as in 
equation ([8]), and rearranging the terms, one finds 



7(r) - 2/3(r) = r 



da? 
dr 



d*T 
dr 



> 



(10) 



as an equivalent, alternative formulation of the GDSAI. Un- 
fortunately, despite the deceptively simple form, for a given 
family of consistent models proving the necessity of the GD- 
SAI from equation (|10[) is not easier than working directly 
on phase-space. 



3 RESULTS 

In this Section we apply Criterion 1 to two families of mod- 
els, and we show that also these systems not belonging to 
the generalized Cuddeford family obey the GDSAI in case 
of phase-space consistency. However, in order to illustrate 
the power of the new result, with the aid of Criterion 1 we 
first re-derive almost immediately the results obtained in 



CMlOa with some lengthy algebra for the generalized Cud- 
deford models. 



3.1 A new proof of the GDSAI for generalized 
Cuddeford systems 

As described in CMlOa, the DF of each component of gen- 
eralized Cuddeford systems is described by the sum of an 
arbitrary number of Cuddeford (1991) DFs with arbitrary 
positive weights Wi, 



f=J 



2Q E 



wMQi), Qi = £ 



2rV 



(11) 



and possibly different anisotropy radii r al , but same h func- 
tion and angular momentum exponent a. Here, a is a real 
number > — 1, and h(Qi) — for Qi < 0. Note that 
multi-component OM models (CP92, CM09), and constant 
anisotropy models are special cases of equation (|11[) . ob- 
tained for i=l and a = 0, and i = 1 and r a = oo, respec- 
tively. 

As shown in CMlOa, the spatial density associated with 
the DF (|ll|l is a factorized function as in equation (poj), where 



A(r) 



(27r) 3/2 r(q+l) y 

V ' fa + 32 ^ 1 



(a + 3 /2) ^ (i + rV^+n 



(* T - Q) 



Q + l/2 



h(Q)dQ, 



(12) 



(13) 



and r(i) is the complete gamma function. Therefore, the 
anisotropy parameter /3, derived in CMlOa by direct com- 
putation, can now be obtained from equations (|12|) and l[7|. 
and Criterion 1 can be applied to the family of generalized 
Cuddeford models. In CMlOa it was shown with some alge- 
bra (see equations [19] and [25] therein) that for a > —1/2 
the first of the necessary conditions for phase-space consis- 
tency can be rewritten as the GDSAI. A more elegant proof 
can now be obtained by using Criterion 1. In fact, direct dif- 
ferentiation of equation (23) shows that B is a monotonically 
increasing function of V&t whenever h > and a > —1/2. An 
interesting situation arises for a = —1/2. In this case equa- 
tion (28) in CMlOa shows that the GDSAI is both necessary 
and sufficient condition for phase-space consistency: in the 
present context dB/d^T = ft($T), and so the monotonicity 
of B is ensured if and only if h is positive, confirming the 
CMlOa result. We are left with the case -1 < a < -1/2. 
As shown in CMlOa, in this case only a sufficient condition, 
coincident with the GDSAI, is available. What are the in- 
formation that can be derived from Criterion 1 in this case? 
Indeed, for — 1 < a < —1/2, it is not possible to compute 
the derivative of equation (23) directly, and so to test the 
monotonicity of B. However, this problem can be circum- 
vented by first integrating by parts, and then performing 
differentiation, obtaining 



dB 
d* T 



= *° +1/2 /i(0) + 



(*t -Q) a+1/2 h'(Q)dQ. (14) 



From the expression above, we conclude that the GDSAI is 
again necessary for phase-space consistency in all general- 
ized multi-component Cuddeford systems with — 1 < a < 
-1/2 and h'(Q) > 0. 



4 Ciotti & Morganti 



Summarizing, for a > — 1/2 the GDSAI is necessary for 
consistency, for a = —1/2 it is equivalent (i.e., necessary 
and sufficient), and for —1 < a < —1/2 it is just suffi- 
cient, but also necessary if h'(Q) > 0. Of course, whenever 
the last condition is not satisfied, the possibility to build 
a counter-example to the conjecture that the GDSAI is a 
universal necessary condition for consistency remains open. 
If one could prove that for — l<a<— 1/2 the inequality is 
necessary whenever h > 0, then one would conclude that in 
multi-component generalized Cuddeford models the GDSAI 
is equivalent to consistency for — 1 < a < —1/2. 



3.2 The Baes & van Hese (2007) anisotropic 
models 

By using Criterion 1 we now show that the GDSAI is obeyed 
by another (quite large) family of model. 

Baes & van Hese (2007, Section 4.2) considered the fam- 
ily of augmented density profiles 



P = Po 



(3 



-2/30 



r 2<5 

' a 



\J/ S 

= 1 ,M> 



where 8 > 0, q < 0, s > 0, < * < *o, and * is the 
value of the central relative potential. From equation (0) 
the anisotropy profile is 



A) + /Mr/r- a ) 2 
1 + (r/r a ) M 



(16) 



so that Po and fioa are the values of the orbital anisotropy 
parameter at small and large radii, respectively; since the 
anisotropy parameter /3 cannot exceed the value of 1 in case 
of positive DF, both /3q and /3oo are < 1. 

Note that the requirement of finite total mass for the 
density distribution (|15[) . i.e. ty(r) ~ 1/r for r — > oo, trans- 
lates into the condition p+2p x > 3, and from the limitation 
on /3oo it follows that p > 1. Now, it is easy to show that 
dB(*)/d* > for < * < * , when s > 0, q < 0, and 
p > 1, and thus Criterion 1 ensures that the GDSAI is nec- 
essarily obeyed also by these profiles. In the particular case 
q = 0, it is possible to prove again the GDSAI by using 
equation (|10p and considering the expression of the radial 
velocity dispersion given in equation (37) of Baes & van 
Hese (2007). Note also that the analytical DF associated to 
these systems is 



k 



(17) 



where the gu are hypergeometric functions (see equa- 
tions [41] and [42] in Baes & van Hese 2007), and so in 
this case the DF is not of the Cuddeford generalized family. 

We finally note that the two models discussed analyti- 
cally in CMlOb, i.e. the Dejonghe (1987) anisotropic Plum- 
mer model, and the Baes & Dejonghe (2002) anisotropic 
Hernquist model, are both special cases of the family (fT5|) . 
and therefore the fact that they also obey the GDSAI is just 
a special case of the result of this Section. 



3.3 The Cuddeford &: Louis (1995) anisotropic 
polytropes 

We finally consider the family of models introduced by Cud- 
deford & Louis (1995). At variance with the previous mod- 
els, these systems are not introduced by using the augmented 
density technique, but from their DF, in a way similar to the 
models in Section 3.1, so that Criterion 1 cannot be applied 
directly. Their DF is 



f(£,J)=£ q -'h(k) 



J 



2ri£' 



(18) 



where £ > 0, and r a is the anisotropy radius. These models 
are known as anisotropic polytropes, and the special case 
h(k) = (1 + k) a was studied by Louis (1993). The formulae 
for p and p r were already obtained also for the more general 
case of / = g(£)h(k), and here we just report the result for 
the case (Tig)) : 



p = 2 3/2 7rB(g,l/2)*«~ i/ ^ 



g-i/2 ,-i f°° h(k)dk 



o (* + «7)«' 



Pl = 2 5/2 7rB(g,3/2)*^ +i/Z 77 



q+i/2 q -i f°° h(k)dk 



o (k + v) q ' 



(19) 
(20) 



where B is the complete Beta function and n — r 2 /r 2 . In- 



cidentally, the validity of equations ()19[) and (20) can be 
easily checked by using equation Q. Note that the conver- 
gence of the energy integral in equation (|19[) requires q > 
near £ = 0. As p r and p axe in factorized form (see also 
Dejonghe 1986, Sect. 1.7.3), we can apply Criterion 1: since 
p oc ^ -1//2 , the GDSAI is satisfied whenever q > 1/2. This 



result can be obtained immediately also from equation (I10[) . 
as the radial velocity dispersion has the remarkably simple 
expression 



2 Pr 

ov = — 



2* T 

2g + l' 



(21) 



The situation is less straightforward in the interval < q < 
1/2. In fact, if a consistent model exists in this interval, then 
it will represent a case of violation of the GDSAI. Therefore, 
it is natural to ask whether it is possible to construct a con- 
sistent dynamical model with q > 0, but violating the GD- 
SAI (q < 1/2). We are not able to answer this question in 
general, however we note that in the self-consistent case of 
finite total mass in which $t ~ 1/r for r — s> oo, volume inte- 
gration of equation (|19p and successive inversion of order of 
integration shows that q > 3/2 is necessary in order to have 
a finite mass for the component under scrutiny (with possi- 
ble further restrictions dependent on the asympotic nature 
of the h function). In other words, there are not consistent 
anisotropic polytropes with finite mass and q < 3/2. There- 
fore the GDSAI is a necessary condition for anisotropic poly- 
tropes in generic external potential when q > 1/2, and for 
all finite mass self-consistent models. It remains open the 
possibility of existence of consistent anisotropic polytropes 
(of infinite mass) with < q < 1/2, which would violate the 
GDSAI. 

Actually the discussion above can be generalized, sim- 
ilarly to what we did in CMlOa when we constructed the 
family of generalized Cuddeford anisotropic systems. In fact, 
we now consider the family of generalized anisotropic poly- 
tropes with DF 
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f(£,J) = £ q ~ 2 Y, w Mk,), k * = ^T£> ( 22 ) 

with different anisotropy radii r a i and positive weights Wi. 
With the definition rji — r 2 /r ai , equations (|19f) and (20) 
become 

p = 2^\B(q, 1/2)^ £ ^ f (g^> (23) 

Pl = 2^B(q, 3/2)*^^ ,-i (24) 

Jo (& + W 9 

and it is immediate to see that all the previous conclusions 
hold also for these general models, even though their DF is 
not of the family JT5J. 

4 DISCUSSION AND CONCLUSIONS 

After the discovery of the "central cusp-anisotropy theo- 
rem" (AE06), an inequality relating the central value of the 
density slope and the anisotropy parameter in consistent 
stellar systems, successive investigations (CM09, CMlOa) 
unexpectedly revealed that the density slope-anisotropy in- 
equality holds not only at the center, but at all radii in a 
very large class of spherical systems (the generalized multi- 
component Cuddeford systems), whenever the phase-space 
distribution function is positive. We call this latter inequal- 
ity the Global Density Slope Anisotropy Inequality (GDSAI). 

In absence of known counter-examples, i.e. two- 
integrals stellar systems with positive DF but violating the 
GDSAI, in this paper we focused on the possibility that the 
GDSAI is actually universal, i.e. it is necessarily obeyed by 
all spherically symmetric, two-integrals systems with pos- 
itive DF. If such conjecture proved true, it would be re- 
markable not only from a theoretical point of view, but also 
for applications. In fact, it could be used to reduce the im- 
pact of mass-anisotropy degeneracy in observational works, 
as orbital anisotropy would be in some sense controlled by 
the local density slope of the stellar distribution in galaxies 
(in the inner regions where 7 < 2). While a proof of this 
conjecture is still missing, some relevant advance has been 
made: in particular, we obtained a new criterion that allows 
us not only to prove in a simple, very elegant, and unified 
way all the previously known results, but also to investigate 
new families of multi-component models. The main results 
of this paper can be summarized as follows: 

(i) By using two previously known and fully general iden- 
tities relating the density profile, the anisotropy profile, and 
the radial pressure in two-integrals systems, specialized to 
the case of factorized systems, we found a very simple con- 
dition equivalent to the GDSAI, namely that the potential 
dependent function in the augmented density is monotoni- 
cally increasing. 

(ii) As a first application of the new condition, we showed 
that all the previous cases, each of them proved with "ad- 
hoc" analysis, are in fact all simple cases of the new relation, 
and the proof is almost immediate. 

(iii) The new criterion is then applied to extend the va- 
lidity of the GDSAI to other models, namely the Baes & van 
Hese (2007) anisotropic models, and the Cuddeford & Louis 
(1995) anisotropic polytropes. This latter family is extended 



to generalized multi-component anisotropic polytropes, and 
it is shown that also in this case the GDSAI holds. 

(iv) As we are not able to show that the monotonicity 
of the potential dependent function in the factorized aug- 
mented density is necessarily monotonic whenever the DF 
is positive, our investigation leaves open the possibility that 
some spherical systems with positive phase-space distribu- 
tion function may violate the GDSAI. Such examples, if they 
exist, may be found (in the class of models studied so far) 
only in the family of Cuddeford models with angular mo- 
mentum exponent in the range —1 < a < —1/2, or in the 
family of Cuddeford & Louis (1995) anisotropic polytropes 
with infinite mass and < q < 1/2. 

Finally, we conclude by noticing that the proof of the 
general validity of the GDSAI could be obtained by using 
equations Q and (5) in all their generality. At the present 
stage we do not have such a proof; however, neither counter- 
exemples are known to us. Moreover, we note that sup- 
porting arguments are provided by numerical simulations 
of N-body systems, whose end-products show correlations 
between f3 and 7 (e.g., see Hansen & Moore 2006, Mamon 
et al. 2006). In addition, Michele Trenti kindly provided us 
with a large set of numerically computed /„ models (Bertin 
& Trenti 2003) and all of them, without exception, satisfy 
the GDSAI. We also verified numerically that the GDSAI is 
satisfied by a large set of radially anisotropic Hernquist and 
Jaffe models with quasi-separable DF constructed by Ger- 
hard (1991) by using the so-called h a circularity function. 
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